Universal scaling relations in molecular superconductors 
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Scaling relations between the superconducting transition temperature T c , the superfluid stiffness 
p s and the normal state conductivity ao(T c ) are identified within the class of molecular superconduc- 
tors. These new scaling properties hold as T c varies over two orders of magnitude for materials with 
differing dimensionality and contrasting molecular structure, and are dramatically different from the 
equivalent scaling properties observed within the family of cuprate superconductors. These scaling 
relations place strong constraints on theories for molecular superconductivity. 

PACS numbers: 76.75.+i, 74.25.Nf, 74.25.Fy, 74.70.Kn, 74.20.De 
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Understanding the phenomenon of superconductivity, 
now observed in quite disparate systems, such as metallic 
elements, cuprates and molecular metals, involves search- 
ing for universal trends across different materials, which 
might provide pointers towards the underlying mecha- 
nisms. One such trend is the linear scaling between the 
superconducting transition temperature (T c ) and the su- 
perfluid stiffness (p s = c 2 /A 2 , where A is the London 
penetration depth), first identified by Uemura et al for 
the underdoped cuprates Q. Recently, scaling relations 
between p s and the normal state conductivity a have 
also been suggested and a linear relation between p s 
and the product <7q(T c )T c was demonstrated for cuprates 
and some elemental superconductors |2j. Here we show 
that these specific linear scaling relations do not hold for 
molecular superconductors, but a different form of power- 
law scaling is found to link p s , <Jo(T c ) and T c . These 
scaling properties hold as T c varies over several orders 
of magnitude for materials with differing dimensionality 
and contrasting molecular structure, and the scaling is 
dramatically different from that of the cuprates. Our 
findings have considerable implications for the theory of 
superconductivity in molecular systems. 

Molecular superconductors are generally regarded as 
members of the wider group of 'exotic' superconductors 
that have attracted much research effort in recent years. 
However, the number of different examples of molecu- 
lar superconductors is now sufficiently large that system- 
atic studies of their properties may be made indepen- 
dently of the other non-molecular exotic superconduc- 
tors. A general feature of all these exotic superconduc- 
tors is the large carrier scattering rate observed in the 
normal state |3J leading to a picture of them as 'bad 
metals' 0. The scattering rate at temperatures near 
T c may have particular relevance for the superconductiv- 
ity, since it is expected that similar carrier interaction 
mechanisms would be dominant in the normal state re- 
sistance and in the pairing of carriers that leads to the 
formation of the superconducting state. It is therefore 
useful to study the correlation between o"o(T c ) and su- 
perconducting parameters such as T c and p s . Fig. 1(a) 



shows p s /c 2 (= ne 2 /m*e c 2 ) and Tc derived from /zSR 
measurements in the vortex state |f| , plotted against 
<7o(T c ) in the highest conductivity direction for a series of 
molecular superconductors. The materials range from a 
highly anisotropic quasi-one-dimcnsional (qlD) organic 
superconductor ((TMTSF^CIO^, through systems of 
two-dimensional (2D) layered organic superconductors 
(BETS and ET salts) to examples of three-dimensional 
(3D) fulleride superconductors; full details are listed in 

Table 1 0i @M m E H M 13 13 13 13 13 IH • 

Note that the parameter values vary over several orders 
of magnitude, which is important for successful deter- 
mination of scaling properties. We find that p s and T c 
are related to <7q(T c ) by power laws of the form cr™ with 
m = -1.05(3) for T c and m = -0.77(3) for p s ; in both 
cases there is a decrease in the strength of the supercon- 
ducting property with increasing conductivity. For com- 
parison, Fig. 1(b) shows a similar plot using the data on 
cuprates and the elements Nb and Pb recently reported 
by Homes et al 0] . Here the overall trend for T c is much 
less clear and the trend for p s shows a broad increase, 
opposite to that of the molecular superconductors, with 
the positive exponent m ~ 0.75. This difference in the 
p s -(Jo scaling between cuprates and molecular systems 
in the high conductivity direction contrasts with the re- 
ported similarity in scaling behaviour between cuprates 
and organics in the low conductivity intcrplane direction 
|2l| . corresponding to m = 0.85 on a p s -<Jo plot. This 
value is similar to the m ~ 0.75 seen for the cuprate high 
conductivity direction in Fig. 1(b). We thus have a situ- 
ation in which the high and low conductivity directions 
in cuprates, along with the low conductivity direction 
in layered organics, all share a similar scaling property 
where p s increases with increasing oq, whereas for the 
high conductivity direction in the molecular supercon- 
ductors p s behaves quite differently in decreasing with 
increasing co- 
in the case of the molecular systems, the different 
power laws seen for p s and T c against cro in Fig. 1(a) im- 
ply that the scaling between them will not be of the linear 
Uemura formQ but will follow another power law. Fig. 2 
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Label Material 


T c 


A 


(To(Tc) 






(K) 


w 


(10 3 Scm- 1 ) 


1 


K-BETS 2 GaCl 4 


0.16(2) [JJ 


2.3(1) [7] 


250(25) [13J 


2 


(TMTSF) 2 C10 4 


1.1(1) [8j 


1.27(3) [8J 


39(6) (a-axis) [14] 


3 


a-ET 2 NH 4 Hg(SCN) 4 


1.1(1) [9] 


1.1(1) [9] 


36(6) [15] 


4 


/?-ET 2 IBr 2 


2.2(1) [9] 


0.90(3) [9J 


26(2) [16] 


5 


A-BETS 2 GaCl 4 


5.5(1) [7] 


0.72(2) [7] 


11(1) [17] 


6 


k-ET 2 Cu(NCS) 2 


9-2(2) [9,^0] 0.54(2) [9, 6(1) [18] 


7 


K3C60 


18.9(1) yjj 


0.48(2) [11] 


2.9(9) [19] 


8 


Rb 3 C 6 o 


29.3(1) [12] 


0.42(2) [12] 


2.5(6) [20] 



TABLE I: Parameter values for the molecular superconductors. Values for T c and A are derived simultaneously from muon 
spin rotation studies in the vortex state. A corresponds to the estimated zero temperature value A(0). o~o(T c ) is the normal 
state conductivity in the most highly conducting direction. The conductivity is derived from reported multi-contact resistance 
measurements in the case of the organics, single-domain STM measurements for K3C60 and far-infrared reflectivity for RbaCeo- 
Estimated uncertainties in the least significant digit are shown in brackets after each value. 
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FIG. 1: (a) The inverse square of the penetration depth 
1/ A 2 (filled circles, left hand scale) and T c (filled triangles, 
right hand scale) plotted against ao(T c ) the normal state con- 
ductivity just above T c in the most highly conducting direc- 
tion. The key to the molecular superconductors is listed in 
Table 1. (b) The data of Homes et al on cuprates and 
elements for comparison; 1/ A 2 (open circles, left hand scale) 
and T c (open triangles, right hand scale). From this data it 
can be seen that 1/ A 2 exhibits a power law dependence on oq 
that is completely opposite to that of the molecular systems. 



shows the Uemura plot of T c against 1/ A 2 in log-log form 
where it can be seen that T c follows p" L with the fitted 
value m = 1.44(3). This approximate scaling of T c with 
pV" 2 1 or equivalently A~ 3 , in 2D organic superconduc- 



tors was noted previously and discussed in terms of the 
2D physics of layered superconductors 0, HH • How- 
ever, it now appears that the scaling relations between 
T c , p s and do are more universal, encompassing examples 
of qlD and 3D molecular superconductors alongside the 
2D systems. The non-linear scaling between T c and p s 
in the molecular case is much harder to understand than 
the linear scaling seen in the cuprates. In the cuprates 
the carrier density n is directly controlled by the doping 
level; in the underdoped regime p s is directly propor- 
tional to n and T c has been suggested to be linked to p s 
either through Bose-Einstein condensation of preformed 
pairs [24| or through a mechanism in which phase fluctu- 
ations of the superconducting order parameter determine 
T c [25j. In contrast, for the molecular systems n is fixed 
by the unit cell size and stoichiometry of the crystal struc- 
ture and varies only little across the range of materials, 
whose superconducting parameters are nevertheless vary- 
ing across several orders of magnitude. Differences in the 
superconducting properties must then arise entirely from 
differences in the details of the electronic many-body in- 
teractions. 

Further evidence for fundamentally different behaviour 
between molecular and non-molecular superconductors is 
seen when an attempt is made to search for linear scaling 
between p s and the product (To(T c )T c , of the form that 
was recently demonstrated by Homes et al 0- Fig- 3(a) 
shows that such a simple linear scaling does not occur 
for the molecular superconductors. The linear behaviour 
seen for the non-molecular systems can be understood 
from applying the Ferrell-Glover-Tinkham sum rule for 
the real part of the frequency dependent conductivity 
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where <j(uj) takes the Drude form <7 /(l + (uj/T) 2 ), with 
(Tq = ne 2 /(m*T) and T being the scattering rate. In 
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FIG. 2: Log-log Uemura plot of T c against 1/A 2 . Data for the 
cuprates and elemental superconductors tabulated by Homes 
et al are also shown for comparison. For the molecular 
systems a scaling close to f?J 2 is observed, rather than the 
linear p s scaling seen for the cuprates. 



the case where T is significantly smaller than the fre- 
quency corresponding to the superconducting energy gap 
2A/7i, the whole free carrier spectrum is redistributed to 
zero frequency to give the superfluid response peak, i.e. 
eojOs = fo(7c)r. If, on the other hand, 2A/h is signif- 
icantly smaller than T, then the normal state conduc- 
tivity is independent of frequency in the gap region, i.e. 
a(u>,T c ) = <7rj(T c ); in this case, as the superconducting 
gap forms, an area of the conductivity spectrum with 
frequency width 2A/h and height <7o is redistributed to 
zero frequency to give the superfluid response peak. This 
leads to the following expression for eop s - 
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where r/ — 2A/fesT c . In Fig. 3(a) the dashed line 
shows Eqn.[3plotted taking the weak-coupling BCS limit 
r\ = 3.53 as a reference; this is seen to describe the gen- 
eral behaviour of the non-molecular data quite well. The 
effective value of ry derived from the data using Eqn.[2is 
shown in Fig. 3(b), which reveals the considerable vari- 
ation among the molecular systems. Note that Eqn. [3 
was derived assuming that the ratio of carrier density to 
effective mass is the same in the normal and supercon- 
ducting states. If, however, this assumption is relaxed 
then the effective gap ratio observed in this plot becomes 
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where the subscripts s and n refer to the superconduct- 
ing and normal states respectively. Strong coupling can 
increase r] over the BCS value via the first term of Eqn.01 
however if Eqn|21 is applicable to the molecular systems 
then the reduced values of r\ seen for many cases would 



require a contribution from at least one of the other two 
terms, i.e. the superconducting carrier density would 
need to be less than normal state carrier density or the 
effective mass of superconducting carriers would have to 
be larger than that of the carriers in the normal state. 
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FIG. 3: (a) Plot of 1/A 2 against the product of T c and a (T c ), 
following Homes et al Q. For the molecular superconductors 
the data collapse onto a narrow ordinate range due to the 
inverse scaling between T c and <to demonstrated in Fig. 1(a). 
Open circles show the data of Homes et al along with 
the linear fit (solid line). The dashed line shows the scal- 
ing expected for a weak coupling BCS superconductor in the 
high scattering rate limit (Eqn|5J. (b) The data expressed as 
an effective gap parameter r\ = (2A/kBT c )(n a /n n )(m^/m* ). 
Whereas the cuprates and elements are grouped around a 
value of rj just above the BCS limit (dashed line) and com- 
parable to the gap ratios seen using other techniques, the 
molecular systems cover a wide range of rj values, both above 
and below the BCS limit. 

Another way to look at the data is to plot the ratio 
en/Os/cn which gives a measure of the effective frequency 
width T c of the normal state conductivity spectrum that 
provides the superfluid response (Fig. 4). This will be de- 
termined either by T itself or by (2/7r)(2A/7i), whichever 
is the smaller. For the non-molecular systems the effec- 
tive width follows the linear T c dependence expected if 
it is proportional to either a BCS-type gap or a T-linear 
scattering rate. In contrast, for the molecular systems 
r c follows a steeper power law T" with the fitted value 
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a = 1.58(5). This behaviour suggests that the molec- 
ular systems are in the low-scattering-rate limit where 
r < (2/n)2A/h and T c follows V. We note that the fit- 
ted power law for T c is also broadly consistent with the 
temperature dependence of the scattering rate deduced 
from the temperature dependent resistance of individ- 
ual examples of the molecular metals; measurements for 
molecular metals just above T c generally show power-law 
exponents in the region 1.5 to 2 El El El El El El El- 
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FIG. 4: Plot of p s normalised by cro(T c )/eo to give the fre- 
quency F . This reflects the effective width of the normal 
state conductivity spectrum, <j(ui), that has condensed into 
the superfluid peak at ui — 0. For the cuprates and elements, 
r e appears to follow a linear temperature dependence. In 
contrast, for the molecular superconductors the steeper tem- 
perature dependence T c oc j^ 08 ' 5 ) [ s se en (solid line). r e for 
the cuprates and elements and higher T c molecular supercon- 
ductors is seen to lie between the BCS weak-coupling limit 
(dashed line, Eqn. |2Jl and the Planckian scattering rate limit 
3] (dotted line). 

The scaling relations for the molecular superconduc- 
tors highlighted here suggest that there are features of 
their electronic properties that are common, despite the 
various materials having quite different dimensionality 
and Fermi surface topology. The simplicity of the scaling 
also suggests that it is being controlled by a single domi- 
nant parameter, such as the ratio of the electron correla- 
tion energy on a molecule U to the electronic bandwidth 
W. U/W also controls the proximity of the Mott insu- 
lator (MI) phase; being close to the MI phase supports 
higher T c (as seen in studies of the K-phase ET salts un- 
der pressure [13, HI) and the inverse relation between 
T c and <7o follows naturally from this. However stan- 
dard approaches to modelling the crossover between MI 
and superconducting phases predict that the enhanced T c 



near the MI phase is accompanied by a depressed p s [29j , 
exactly the opposite of what is observed experimentally. 
Identification of new theoretical models that match the 
observed scaling behaviour is clearly necessary; finding 
such models should lead to significant progress in under- 
standing superconductivity in molecular systems. 
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